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We present our calculations of the electromagnetic form factor of pions. We explore the prop- 
erties of pion form factor at momentum transfer larger than previous studies by including more 
combinations of source and sink momenta and using more configurations. We fit our results using 
vector meson dominance (VMD) hypothesis. 
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1. Motivation 

The pion form factor is often considered a good observable for studying the onset of the pertur- 
bative QCD (pQCD)regime in exclusive processes. There are several reasons: First, the asymptotic 
forms of the pion form factor at both large and small are known. At large it scales as 
[1,2,3,4,5,6, 7,8,9] 



Q2 



(1.1) 



while at small Q^, the pion form factor can be well described by the Vector Meson Dominance 
(VMD) Model [10, 11, 12] 



1 



for < m^MD 



(1.2) 



l + eV'"vMD 

Therefore at some there must be a transition from the VMD behavior to the large scaling 
predicted by pQCD. Since the pion is the lightest hadron, the transition is expected to occur at 
lower than heavier hadrons, which makes it relatively easier to probe by both experiments and 
Lattice QCD (LQCD). Finally, there is no disconnected diagram on the lattice for the pion form 
factor. Thus the calculation is pretty straightforward. Previous LQCD studies on the pion form 
factor can be found in [13, 14, 15, 16] and the references therein. 

The current results from various experiments are shown in Fig. 1, including the latest results 
from Jefferson Lab (JLab) experiments E93-021 [17, 18] and EOl-004 [19]. As we can indicate 
from the figure, the data points around ~ 2GeV^ start to show some hints of deviation from the 
VMD fit. This is the energy regime we would like to explore in our study. 
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Figure 1: Summary of experimental data for the pion electromagnetic form factor. The two points with 
open circles are the latest data from the Jefferson Lab (JLab). Shaded regions are expected sensitivities of 
future experiments. 
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2. Lattice Techniques 

In this section we explain the techniques we used in our lattice calculations, namely the sequen- 
tial source method (for calculating the quark propagator) and the ratio method (for the correlation 
functions). The pion electromagnetic form factor is obtained in LQCD by placing a pion creation 
operator (the "source") at Euclidean time f,- with momentum p,, a pion annihilation operator (the 
"sink") at tf with momentum pf, and a current insertion at time t with momentum transfer q, as 
shown in Fig. 2. The standard quark propagator calculation provides the two propagator lines that 
originate from f,-, the remaining propagator from is obtained by the sequential source method: 
completely specify the quantum numbers and p f at the sink, and contract the propagator from f, to 
tf with the annihilation operator to serve as the source vector of a second, sequential propagator 
inversion. The advantage of using the sequential source method is that various currents with differ- 
ent can be inserted at time t without additional matrix inversions. The largest available lies 
in Breit frame (pf = —pi). 




Figure 2: The quark propagators used to compute the pion form factor. 



To obtain a simple expression on the lattice, we construct the pion form factor using the ratio 
method. The pion form factor F{Q^) is defined as 



continuum 

= Zv {K{pf) \V^{0)\K{pi)) = F{Q^){pi + pf)^ 



(2.1) 



where V^{x) is the chosen vector current. We can extract F{Q^) from some ratio of the three-point 
correlation function and the two-point functions. The three-point function can be written as 



^ilK{h,t,tf,pi,pf)=a'^ £ e 



-i(xf-x)-pf 



xe 



(2.2) 
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where 0's are operators with pion quantum numbers; A € {L,S) and B G {L,S) denote either 
"local"(Z.) or "smeared" (5). Inserting complete sets of hadron states and requiring tj <^ t <^ tf, 
gives 



^ifiAti,t,tf,Pi,Pf) ^ {0\(()B{x)\7l{pf)) 
X {K{pf) \V^{x)\K{pi)) Uipi) 0](x) 



4E„{pf)Ej,{pi) 
Similarly for the two-point correlator, 



,-{tf-t)E„{pf) -(t-ti)E^{pi) 



(2.3) 



Tii{tutf,p)^mB{^i)\Tl{p)) 



x{k{p) 



IE 



We can obtain F{Q^) from the following ratio 



(2.4) 



^En{pi) +EnPf^ 

where the indices A, B and C can be either L (local) or S (smeared). 



Ti^^{ti,t,Pi)TCB{ti,tf,pf) 

2ZvEn{pf) 



(2.5) 



3. Simulation Details and Results 

We use lattices generated by MILC [20], with volume 20^ x 32 and lattice spacing a = 0.125 
fm. The sea quark mass nisea and the valence quark mass nivai are tuned so that we get the same 
lightest pion mass mji{msea) = iiin^mvai) [21]. The pion operators are fixed at time = 10 and 
tf = 20, and the number of configurations used in this study is 201. We use five different sets of 
sink momenta: ^/ = (0,0,0), (1,0,0), (1, 1,0), (1, 1, 1), and (2,0,0). 

We present our results in terms of the square of the pion charge radius, obtained by the VMD 

fit: 



2 



(3.1) 



as shown in Fig. 3. The first point on the left is from the data set with only zero sink momentum 
(pf = (0,0,0)), and for the next point we combined the data from both pf = (0,0,0) and pf = 
(1,0,0), and for the third point we added in pf = (1,1,0), and so on. 

We can see from Fig. 3 that the error bai^s of r| increases as higher sink momenta ai^e included. 
Since the pion charge radius is related to the slope of F{Q^) at low Q^, we derive (r^) from the data 
set of zero sink momentum pf = (0,0,0) and < IGeV^, and check the consistency between the 
VMD fit and the data above IGeV^ to see if there is any deviation from the VMD model. 
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Figure 3: Pion Form Factor VMD fit for = (0,0,0) to (2,0,0). 

The result of this "consistency check" is presented in Fig. 4, where we plot Q^F{Q^) against 
Q^. While the quantity Q^F{Q^) should approach a constant as predicted by VMD, we can see 
that there are some hints of deviation from the VMD model for points with > 2GeV^. To 
further emphasize this observation, we define ^Q^F{Q^) = Q^F{Q^)icntke — Q^F{Q^)vmd, and 
plot AQ^F{Q^) against in Fig. 5. 
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Figure 4: Consistency between data and the VMD fit from low with zero sink momentum. The three red 
lines represent the VMD fit and its error bars, and the triangle points correspond to the experimental data 
from JLab. 

We also compared this VMD fits from low with that from the Breit frame (wherep/ = —pd, 
for the data points of the Breit frame have relatively small error bars at high Q^. The result is shown 
in Fig. 6. The figure implies that the VMD fit from the Breit frame (the purple one) is about la 
away from the fit of a single zero sink momentum (the blue one), hence exploring further in the 
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Figure 5: ^.Q^F{Q^), as defined in the text. 
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Figure 6: VMD fits with error bands from both data of low Q with zero pj, and from data points in the 
Breit frame. 



Breit frame may be the correct direction for studying the pion form factor at higher momentum 
transfer. 

4. Summary and Outlook 

In this study we have acquired enough lattice data for < IGeV^ to extract a reliable pion 
charge radius r,^. By comparing the VMD fit from data with low and high Q^, we've started 
to see some hints of discrepancy between data points at different momentum transfer, which may 
indicate the transition from the VMD model to pQCD, the goal we are seeking for. We also found 
that the VMD fit from the Breit frame is about 1 a away from the fit of a single zero sink momentum, 
and we infer that we may explore further in the high regime by studying the data from the Breit 
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frame. We are generating four times more data to shrink the error bars in the Breit frame in the 
hope of a clearer and stronger evidence of the transition into the perturbative QCD regime. 

In the meantime, the JLQCD Collaboration has also reported their calculation of the pion 
form factor based on all-to-all propagators. Interested readers may find details in their upcoming 
publication [22]. 
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